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We present the non-linear theory of shock acceleration applied to SNRs ex- 
panding into partially neutral plasma. Using this theory we show how the 
Balmcr lines detected from young SNRs can be used to test the efficiency of 
shocks in the production of cosmic rays. In particular we investigate the effect 
of charge-exchange between protons and neutral hydrogen occurring in the pre- 
cursor formed ahead of the shock. In this precursor the CR pressure accelerate 
the ionized component of the plasma and a relative velocity between protons 
and neutral hydrogen is established. On the other hand the charge-exchange 
process tends to equilibrate ions and neutrals resulting in the heating of both 
components. We show that even when the shock converts only a few per cent 
of the total bulk kinetic energy into CRs, the heating is efficient enough to 
produce a detectable broadening of the narrow Balmer lines emitted by the 
neutral hydrogen. 
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1. Introduction 

Supernova remnants (SNRs) are thought to be the primary sources of 
Galactic cosmic rays (CRs). Diffusive shock acceleration (DSA) occurring 
around the forward shock is considered the most promising mechanism for 
the acceleration of particles up to very high energies. In order to explain 
the flux of CRs observed at the Earth, the shock should convert a frac- 
tion ~ 10 — 30% of the plasm kinetic energy into non thermal particles. 
Indeed, from the theoretical point of view, DSA can be very efficient in 
producing accelerated particles, but a clear confirmation of this prediction 
still lacks, even if many circumstantial evidences have been collected (see 
e.g. Ref. 1). A possible way to investigate the efficiency of CRs produc- 
tion is through the analysis of the Balmer lines associated with shocks of 
SNRs which propagate in a partially ionized plasma. The hydrogen lines 
of a Balmer shock consist of two superimposed components: the narrow 
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component is emitted by neutral hydrogen after entering the shock front 
and the broad component by hot protons after undergoing charge exchange 
with incoming neutral hydrogen atoms. The width of the broad component 
reflects the proton temperature behind the shock front, while the width of 
the narrow component provides a direct measurement of the temperature 
upstream of the shock. 2 ' 3 

We stress that, when the shock is not modified by the presence of CRs, 
the charge exchange process and the Balmer emission both occur only 
downstream of the shock. On the other hand, when CRs are accelerated 
efficiently, the shock structure is modified and a precursor is generated up- 
stream: the CRs pressure slows down and compresses the ionized plasma. 
Even if the neutral component is not directly affected by the CR pressure, a 
relative velocity between ions and neutrals is now established. Hence charge 
exchange can also occur upstream of the shock, resulting in the pre-heating 
of neutral hydrogen. We do expect two remarkable consequences: 1) the 
narrow Balmer lines can be emitted also from the upstream region and 2) 
the typical width of these lines becomes larger with respect to the case 
without CRs, because the hydrogen temperature increases. 

A further consequence of efficient acceleration is that the temperature 
of the shocked ions downstream is lower with respect to the case with no 
acceleration. This occurs because a non negligible fraction of the shock 
kinetic energy is converted into CRs, rather than into thermal energy. As 
a consequence also the width of the broad Balmer lines is affected, being 
reduced with respect to the case with no acceleration. 

Remarkably all these aspects have been observed in some Balmer- 
dominatcd shocks:f ) from a region of the Tycho remnant, the emission of 
narrow Balmer lines has been detected upstream of the shock; 4 2) in several 
SNRs the narrow Balmer lines present a width incompatible with the typical 
temperature of the interstellar medium; 5 3) in two different cases where the 
Balmer emission has been detected, i.e. RCW 86 6 and SNR 0509-67. 5, 7 the 
width of the broad lines led to a downstream temperature which is lower 
than that estimated form the measurement of the shock proper motion, 
suggesting that a fraction of the total energy is converted into CRs. 

In order to quantify all these effects simultaneously, in this paper we 
present the solution for the stationary case of DSA when the shock prop- 
agates into a partially ionized plasma. We use a semi-analytical method, 
similar to those used in Ref. 8, specialized for plane shock geometry. 
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2. Interaction between ions and neutrals in the precursor 

The process of charge exchange can occur in the precursor only if its typical 
length scale is smaller than the precursor length. The length of the precursor 
created by CR pressure is of the order of the diffusion length of the particles 
with maximum energy, i.e.: 

L prec ~ D(E mayL )/uv = 3 • IQ 17 E TcV B~qU% 1 cm , (1) 

where the maximum energy is expressed in TeV, the magnetic field in 
fiG and the shock speed is uq = u 8 10 8 cm/s. We assume that particles 
diffuse with a Bohm-like diffusion coefficient, D{p) = (c/3)rt,(p), where 
tl = pc/eB is the Larmor radius and the magnetic field strength is as- 
sumed constant in the precursor. Let call tin the neutral hydrogen density, 
the length scale for a proton to make charge exchange upstream is: 

Lee = u / (n H <r ce (v r ) v r ) ~ 10 15 (u Q /v r )(n H /lcm^y 1 cm . (2) 

We note that when the relative speed between ions and neutral hydrogen 
is v r < 2000 km/s, the charge exchange cross section is roughly constant. 3 
Comparing L ce and L prec we see that a single proton can charge exchange 
several times in the precursor before crossing the shock. In this work we 
neglect the role of ionization because the ionization cross section is much 
smaller than the charge exchange length when v r < 1000 km/s, but we will 
defer a more extensive treatment of ionization to a future paper. 



3. Hydrodynamics with charge-exchange 

The behaviour of a system of neutral hydrogen and protons interacting via 
charge exchange is described by a kinetic equations for each species: 

dfi 9fi 

~dt + Vx ~dx = v )-^'( a; ' v ) ~ PA x > v )fi( x iV) , (3) 

where i, j — H,p and /j is the distribution function. Eq. (3) is written in 
one spatial dimension because we are considering plane shocks propagating 
along the x direction. As a consequence /j is only a function of the position 
x and of the particle velocity. The charge exchange frequency is defined by: 

Pi(x,v) = Jdw\v- w|<7 ce (|v- w|)/i(a;,w). (4) 

Notice that the charge exchange cross section is only a function of the 
relative speed v r = |vjj— v p |. Now we consider only the stationary solution, 
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hence dfi/dt = 0. From the integration of first, second and third moments 
of Eq. (3) we get the fluid equations for neutral hydrogen: 



^[PHU H ]=0, 

d 

[phu 2 h + P H ] = -q m , 



d_ 

dx 



;Phu 3 h 



7 



7-1 



Phu h 



= -1e 



(5) 
(6) 

(7) 



The term q m and q e are the flux of momentum and energy which are trans- 
ferred from neutral hydrogen to protons, i.e: 



q m = J dv H dv p (v H - v p ) a ce {v r )v r f H (v H )f p (v p ) 
q e = J dv H dv p ^(v% - Vp) <T ce (v r )v r f H (v H )fp(v p ) . 



(8) 
(9) 



In order to get an analytic expression for q m and q e we follow Ref. 9. First of 
all we assume that both hydrogen and protons distributions are Maxwellian, 
with a bulk velocity Ui(x), along the x direction, and a temperature Ti(x), 
i.e. fi(x,v) — pi (it u§,.) _3 / 2 e~^ v_u ^ , where the thermal speed is de- 
fined as = \JlPij pi. Moreover we note that for v r < 2000 km/s, the 
charge exchange cross section decreases only logarithmically with increasing 
v r , hence we can approximate a ce as a constant evaluated at the average 
relative speed and pull it out of the integrals. Under this assumptions mo- 
mentum and energy transfer can be expressed as follows: 

q m = <Jce(u)ppPfi(uH - u p ) [u + v TH /hi(H,p) + v Tv /hi{H,p)] ; (10) 



qe = (Jce(u)PppH 



K) + - 2 



Vt h h 2 (p, H) - vlh 2 (H, p) 



+ 



+2(u H - Up) 



(11) 



h!(p,H) fti(H.p) 
where the average relative speed u, and the functions are defined as: 

'A 



u = 



: (v T . + v T .) + (ui -UjY 

\ 97T 

4 ( + { Ui - uj) 2 + —v 2 



7T 



4 Q 64 9 / \ 2 

~ V T, + ~n~ V Tj + K - u j) 
IT 97T J 



1/2 
1/2 
1/2 



(12) 
(13) 
(14) 
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The fluid equations for protons are similar to those for hydrogen, but 
include the contribution due to CRs. Following Ref. 8 it is easy to show 
that the flux conservation of mass, momentum and energy are: 



9 r 







d 

— [ Pp u 2 p + P P + P c ] = q, 
dP c 



d_ 

dx 



7 



7 



1 



-PpXLp 



(15) 
(16) 
(17) 



where P c is the CR pressure which is connected to the CR distribution 
function through P c (x) = 47r/3 / p 2 dppvf c (x,p). Equations (5)-(7) and 
(15)-(17) are coupled through q m and q e and can be solved numerically 
to get the pressures and the bulk velocities of both protons and hydrogen 
as implicit functions of P c . In order to close the system of equations we 
need to solve also the transport equation for the CR distribution function, 
f c . All the system can be solved using iterative techniques similar to those 
used in Ref. 8. We will describe the full procedure in a more detailed work. 

It is worth stressing that the technique used allows one to take into 
account other important effects in the shock acceleration theory, like mag- 
netic field amplification and turbulent heating. For the sake of simplicity 
here we neglect such complications in order to isolate the role of charge 
exchange. 



4. Results and Conclusions 

We consider a typical case of shock with the following parameters: shock 
speed Mo = 2000 km/s, total upstream density po = 1cm -3 , neutral frac- 
tion 50%, upstream temperature T = 10 4 K, upstream magnetic field 
strength B\ = \0pG and maximum momentum of accelerated protons fixed 
to Pmax = 10 4 77j p c. The parameter which regulates the injection efficiency 
is taken £; n j = 3. 95. 8 For these values we get L prec /L ce ~ 150, hence we do 
expect an efficient role of charge exchange, while L prec /Li on ~ 10~ 4 hence 
the ionization can be neglected as stated above. In Fig. 1 we show the ve- 
locity and the temperature profiles in the precursor for both protons and 
hydrogen. The dotted line shows the normalized CR pressure, which at the 
sub-shock position reaches the value P c / pqUq = 0.14, hence we are dealing 
with a mildly efficient shock. The upstream temperatures of protons and 
hydrogen reach the values of 5 • 10 5 K and 4.7 • 10 5 K, respectively. Protons 
are slightly wormer than hydrogen because they also suffer the adiabatic 
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compression due to the CR pressure. The narrow Balmer line width corre- 
sponding to the hydrogen temperature at ths subshock is 46 km/s, to be 
compared with 21 km/s resulting in the absence af a CR precursor. The 
downstream temperature of shocked ions is 6.7 ■ 10 7 K, which produces a 
broad Balmer line with a width of 1774 km/s; in the absence of acceler- 
ation the same width would be 2039 km/s. As already pointed out, this 
difference arises because when the acceleration is efficient a non negligible 
fraction of the total shock kinetic energy is channelled into CRs rather than 
into thermal energy. 
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Fig. 1. Profiles of velocities and temperatures of hydrogen and protons in the precursor 
upstream of the shock. The thick-dotted line shows the CR pressure. All quantities are 
normalized as shown in the legend. The typical precursor length is x max = D(p max )/un. 
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